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THE n.OV/ OF A COJffRESSIBLE FLUID PAST 
A CURVED SURFACE 
3y Carl Kaplan 

SUiaiARY 

An iteration method is employed; to obtain the flow of a com- 
pressible fluid past a curved surface. The first a.pproxiTnation, 
v/hich leads to the Prandtl-Giauert rule, is based on the assuiaption 
that the flov/ differs but little from a pure translation. The 
iteration process then consists in improving this first approxi- 
mation in order that it will apply to a flov; differing frora pure 
translatory motion to a greater degree* The method fails v^hen the 
Ife-ch number of the undisturbed stream reaches unity bwt permits a 
transition from subsonic to supersonic conditions without the ap- 
pearance of a compression shock. The limiting value of the 
undisturbed stream Ifech number, defined as that value at which 
potential floAv no longer exists, is indicated by the apparent di- . 
vergenoe of the power series representing the velocity of the fluid 
at the surface of the solid boundary. 

For small }^ch numbers and for thin shapes, the results ob- 
tained by the iteration process agree v/itVi those obtained by the 
Poggi method. For higher values of the stream Mach number less 
than the critical value, numerical calculations are in agreement 
with the results obtained by von Karman by means of the hodograph 
method. For values of the stream Mach numbc-;r higher than- the 
critical value, the iteration process yields some information ^bout 
the region of flow comprised between the critical stream Ifech num- 
ber and the limiting stream Mach number. 

INTRODUCTION 

When a body is held fixed in a compressible 'f3.uid moving at a 
uniform speed less than, but comparable vdth, that of sound, there 
may be a region near the surface v/here .the velocity of the fluid 
relative to the body exceeds the local velocity of so'Jnd. The flow 
in such cases may be perfectly regular with no indication of shock 
waves. Several' such tjrpes of flov/ have been described by Taylor 



2 



(reference 1) and, more recently, by Gortler (reference 2). In con- 
nection with this type of flow it is important to knov/ when the ir- 
rotational motion ceases to be possible. It is certain that ir- 
rotational motion is no longer possible as soon as the !-^ach number 
of the undisturbed stream reaches a definite value, alv/ays less than 
unity, v/hich depends on the shape of the body. In the past com- 
pressibility shock has often been assumed to occur when the maximum 
velocity of the fluid at the surface of the body equals the local 
velocity of soimd; hov/ever, the papers of Taylor and Gortler (ref- 
erences 1 and 2) question vriiether this is correct for the first ap- ^ 
pearance of a shock wave. In addition, a recent paper by von K^rman 
(reference 5) suggests that the envelope of the Ivfeoh lines in the 
supersonic region of flov^ probably introduces the first shock wave 
in the flow. The stream Mach number at which the envelope of the 
liliach lines first appears may thus be identified v/ith a limiting 
value of the Jlach number. 

The present paper treats the flov/ of a compressible fluid past 
a- curved siarface by means of an iteration process based on that of 
Ackeret (reference U) • The boundary was so c?iosen as to conform 
v;ith the requirements of the method; namely, no stagnation points 
and small variation of the local velocity from that in the stream. 
The process, furthermore, permits values of the stream Mach number 
ranging from zero to the neighborhood of unity. The method proves 
to be quite laborious when more than t-wo stages in the iteration are 
deme.nded; but, because of the importance of the problem, it has been 
thought worth while to perform, the .third step. Most af the details 
of the calciolations have been relegated to appendixes in order not 
to disturb the continuity of the main ideas. (The equations in the 
appendixes have been assigned nuiabers prefixed by letters denoting 
the appendix; for extim.ple, equation (A-5) is the third equation in 
appendix A.) 

THE ITSBATION PROCESS 

The fundamental differential aquation governing the flow of a 
compressible fluid is 

where 



X, Y rectangular Cartesian coordinates in plane of profile 
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V fluid velocity oomponentG along X and Y axes 
c local velocity of sound 

The condition for irrotational motion is that 

bY dX 

and leads to a velocity potential p' defined by 



u 

OX j 



V (2) 



S> I 

CY 



If the body is held fixed in a uxuform stream of velocity U, the 
relation between the local velocity of sound c axid the velocity of 

the fluid + ^2 is given for adiabatic processes by 




(3) 



where 



'O 



velocity of sound in undisturbed stream 



Y ratio of specific heats at constant pressure and constant 

volume 

M I*ich number of undisturbed stream (u/cq) 

Vftth the introduction of a characteristic length s. as unit of 
length and the stream velocity U as xmit of velocity, the various 
quantities thus far defined can be rendered nondimensional. Thus 
X, . Y, u, V, and / denote, respectively, the nondimensional 
quantities X/s, y/s, u/h, v/u, and ^/Us v/hile c and c^ 
retain their original meanings. By use of equation, (2), equa- 
tions (1) and (5) then become, respectively, ' / ■ 



and 



7^ 



* 1 



(5) 



Let t denote a characteristic parameter of the shape, such as 
the thickness coefficient: then, the foil owing expansions are as- 
sumed: 



6x 6x dx 6x 



6y 5y dY dY 



(6) 



1/Vhen these expressions for u and v are introduced into equa- 
tion (1-0, together vath the expression for oVCq'" given by equa- 
tion (5), and when the coefficients of the' various powers of t are 
equated to zero, the following differential equations for jzf^^, j^, 
. . • result: 
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d3r 6y- 



(t-1) 



6^ 

6x2 



o ^, 



6^ 
6X 



6_^ 
6x 



6% 



6Df 6y2 



1 



/gj2f 6^, 62^ 64 c^j^, . 6)21' 624N I 
\6x 6y 6x6y. 6y 6x6r 6y 6x6y/ J 

These differential equations may be put into more faniliar 
forms by introducing a next, set of independent variables x and y, 
where 



X » X 



"I 



(10) 



Y 



Thus, for M < 1, equation (7) is transformed into a Laplace equa- 
tion and equations (8) and (9) into Poisson equations. The solu- 
tion of equation (7) yields the well-kno;vn Praiidtl-Glauert rule. 
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whereas the solutions of equations (8) and (9) provide higher ap- 
proximations to the flov/ of a compressible fluid and thus will ap- 
ply for larger departures from an undisturbed uniform flov/.. 

The procedure to be followed in solving equations (7), (8)* e.nd 
(9) is very simple in principle. The first step is to obtain an 
expression for the velocity potential of the incompressible flow 
past the chosen boundary and to express it as a power series in the 
thickness coefficient t. Then the solution for the first approxi- 
mation J2f] to the compressible flow is easily obtained by analogy 
from the coefficient of the first power of t. The second and 
third approximations ^2 ^^'^ ^3 obtained by solving equa- 

tions (8) and (9). The boundary conditions - that the flovj- is 
tangential to the solid boundary and that the disturbance to the 
main flow vanishes at infinity - are satisfied to the same pov/er of 
the thickness coefficient t which is involved in the expression 
for the velocity potential ^* 

Flow past a Curved Surface 

The solid boundary chosen for use in this paoer is a sj/mmetri- 
cal shape ivith cusps at both the leading and the trailing edges, 
thus insuring no stagnation points in a uniform flov/ parallel to 
the axis of symmetry. Appendix A contains the derivation of this 
shape and also the solution for the flow of an incompressible fluid 
past it. Appendixes B, C, and D contain the detailed calculations 
for 02» ^•^'^^ respectively* The final expression for the 

velocity potential p takes the .follovdng form: 

t / -2; -3b , \ 

0 =* cosh £ cos f] + 7-H 3© cos T) - e cos pr]j 

+ t^(Fj cos T) + cos 3'n + P5 cog 5*0) 
+ t-^ Q-^ cos T] + cos J-n + cos 5^1 + cos 7'n 



/ sinh cos rj .^e ^cje 

+ Gq ( = e =*cos IT - e =^ cos 3^1 « e cos 5'Q 

\cosh 2E, - cos 2ri 



-7l ^ 

- 0 cos 7^) 



+ .. . . (11) 
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v^here 

T| elliptic coordinates related to the Cartesian coordinates 
X and y by the equations 



X « cosh § cos Tj 
y « sinh ^ sin r\ 



y (12) 



p « \/l ^ m2 

^1* ^3' ^5 functions of c and of M given by equations (C-18) 

Gi, Gj, G^, Gy, Gq functions of and M given by equations (D-l6), 
(D-17), (D'16), (D-19), and (D-22), respectively 

Equation (11) represents the solution of the fundamental dif- 
ferential equation (1) that satisfies the boundary conditions at the 
surface of the body and at infinity, insofar as the terms inclu- 
sive of the third power of the thickness coefficient t- are con- 
cerned. The coefficients of the various powers of t are exact 
and are valid for all values of the Mach number M from zero up to 
but not including unity. On the other hand, the method of Poggi 
yields the components of the fluid velocity in the form or pov^er 
series in lyT, the coefficients of which are exact and valid for 
the entire range of values of the thickness coefficient t. Ap- 
pendix E contains the solution of the problem of this pap-er by the 
method of Poggi as far as the. terms are concerned. 

Because the iteration process and the Poggi method yield solu- 
tions of the same equation (1) in the form of power series in t 
and in M^, respectively, the two m.ethods must agree in the range 
common to tJiem; that is, the iteration expression for the fluid ve- 
locity at the solid boundary, expanded according to powers of M 
and with all terms containing pov/ers of M higher than the second 
neglected, must agree with the corresponding Poggi result, expanded 
according to powers of t and with all terms containing powers of 
t higher than the third neglected. This calculation is shown in 
detail in appendix F. 

Nvttnerical Applications 

Calculations are now made for the velocity distribution at the 
surface of a bump - that is, a member of the family of shapes derived 
in appendix A - for several values of the Mach number. Because term^s 
involving powers of the thickness coefficient t higher than the 
third have been neglected throughout the present paper, the fluid 
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velocity- q should be expressed in the following form: 

q « 1 + a^t + agt^ + a^t^ + . . • 
From equations (F- 6) it follows easily that 

ai = - cos 2a 
2p 



+ .-i- - -i. + cos 2a 



8p2 16 Up 



+ 3 6 + log - 

16 b2 



cos Ua 



- ^1-p) + 2(3(W + 5Gc + 7G7) 



cos t^a 



cos li-a 



P 



11,(07)^ 



cos 



6a 



> 



(13) 



ilk) 



The expressions for ('^'i)q' ('^)o' ('^)o' (*^)o ^^"'^^ 

the end of appendix P, and table IV shoxvs the calculated values for 
M = 0.50, 0.75. 0.83, and O.S^O. Table V gives the calculated 
values of a-j^, ag, and a^ at various positions along the profile 

for M = 0.50, 0.75, 0.85, and O.9O. With t = 0.10, the expres- 
sions- for the maximum fluid velocity q^^^^^ at the surface can be 
written as follows; 



1 + 0.17321 + 0.0227i». + 0.003i+ii 

1 + 0.25026 + 0. 051+35 + 0.02188 

1 + 0.26895 + o.iiolU + 0.072^0 

1 + 0.3i4l+i2 + 0.19268 + O.U578U 



for M = 0.50 

for M = 0.75 

for M " 0.83 

for M «» 0.90 



(15) 



An examination of the foregoing series shows that qr^yr will probably 
diverge for some value of M in the neighborhood of 'o,83. This 
value of M marks the limit of irrotational potential flow and 
probably indicates the first appearance of a compression shock at 
the solid boundary. Farther on in this section a rule will be for- 
mulated for estimating this limiting value of the stream !/;ach 
number . 



The velocity distribution for a profile of thickness coeffi- 
cient t = 0.10 is calculated by means of table V and equation (13). 
Table VI lists the values of q for M « O.50, O.75, and O.83 and 
figure 1 shows the corresponding graphs. The broken curve repre- 
sents the velocity distribution for t « 0.10 and M * Q.50 calcu- 
lated according to the. Poggi method. (See table HI. ) The curves 
of figure 1 show the agreement betv/een the values of q calculated 
by means of the Poggi and the iteration methods for M ^ O'.^O and 
also the gradual change in cixrvature of the velocity-distribution 
curves in the neighborhood of the leading and the trailing edges as 
the stream Mach number is increased. 

The critical Mach number, defined as that value of the stream 
Mach number at v/hich the local fluid velocity first attains the local 
speed of sound, is calcxilated as follows: - . 

2 ■ 

In equation (5)# {u^ + v^) is put equal to .2-*, or 

^or'~S*r-^) .(16, 

Table VII lists values of q for various values of the Mach num- 
ber M. From equations (ll|) with a = and the expressions for 

(^1)0* (^)o' (^5)0' (^7)0' ^* follows that the maximum velocity 
at the boundary is given by 
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Table II gives the values of A, B, 0, and D obtained from equa- 
tion (D-5) for various values of the Mach number. Table VII lists 
the values of q^j^^x ^'^^ t « 0.10 calculated by means of equa- 
tion (17). Values of q^^^g^^ calculated by the Pogg:i method are 
also given in table VII. Vor low Mach mombers the approximate 
values for q^^^ obtained by the tvvo. methods ar;ree, but for high 
Mach numbers the Poggi method yields values that are too low. 
Figure 2 shows the graphs of q^^ and of q^^^^^ as functions of the 
Mach number. The intersection of the two curves gives the critical 
Mach number. The iteration method (solid line) yields the value 
Mor ^ 0.7li whereas the Poggi method (broken line) yields the value 
Mcr ° 0,788. 

The value of q^j, for iil 0.33 obtained from equation ( l6) 
is 1.1731, The last column in table VIII shows the values of 
q/q ^ for the shape t « 0.10 with M'« O.S?. Values of 
designate a region in which the flow is supersonic whereas values of 
q^/q < 1 characterize the subsonic region. The supersonic region 
of flow is symmetrical with regard to the Y axis and corresponds to 
a calculation given in reference 1 for the flow through a nozzle, in 
which a similar supersonic region of flcm was found at the wall at 
the narrov/est cross section of the tube. 

In order to find .the extent to which the supersonic region 
penetrates the flow, it is sufficient to utilize for the calculation 
only the terms inclusive of t^, since the series for q converge 
rapidly away from the solid boundary. Thus, along the Y axis, the 
follovdng expression for q. is obtained from equations (F-5) and 
(F-U): 
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2p { • cosh I. - 



Then for t « 0.10 and M « 0.85, 

q ^ 1 + 0.i;05l6e"^^+ 0.02753e'^^'^ 

+ -—^ " (-0.02762e"^ + 0.012U2e"^^ - 0.01520e""^^) 

cosh 



The value of ^ for which q = q^^ = 1.17S1 is 0,o8, The corre- 
sponding value of Y obtained from equation (B-l6) is 

Y = — sinh § 
P 

«=0,720 

The supersonic region of flov; thus extends into the fluid a distance 
equal to almost seven times the maximum height of the bump. By use 
of several chosen values of T|, the constant velocity profile 
q z ^ 1«1731 can be plotted by means of equations (F-3) &iid 
(F-U). In figure 1 the lower broken curve represents this profile. 
The region inside the profile is completely supersonic and therefore 
contains real Mach lines. The region outside the profile is every- 
where subsonic and therefore the Mach lines are imaginary. 

The pressure coefficient Cp.^jj is obtained from the expression 

il + i(Y-l)M2(l-Q2) - 1 
^P>M -- , , ^ (^9) 

2 ^ 
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where 

2 0 



and q- is the velocity of the coinpressible fluid, referred to the 
velocity U of the undisturbed streaia. 

Since, throughout this paper, terms involving pov^ers of - t 
higher than the third have been neglected, jj should be expressed 
as a power series in t. Thus, if q from equation (I3) is sub- 
stituted into equation (19)* it follows easily that 



-2a it + [j- (aj^ + 2a2) + ai%2j 



-2^aj^a2+ a^^ + (a^^ + 2a2 ) a^^M* 



2_2-Y _5Ji- 



t" + . . . 



(20) 



With the help of table V, values for Cp j^j along the profile 
can be easily calculated for the case t O.ic and M « 0,83. 
Table IX shows these values of Cp^j/[ .together with corresponding 
values calculated according to the Prandtl-Gls uert rule and the 
von JC^rman method. The Prandtl-Glauert rule is 



^P,] 



%,0 



and the relation obtained by von KaViran (reference 3) is 



C. 



'p,M 



p,0 



^/rri2 



1 + \/r- ^ 



Figure 5 shows the graphs of the various calculated results. It is 
observed that the results of the Prandtl-Glauert and the von .KaVmc'^n 
methods differ considerably from the results of ' the iteration method. 
The "reasons for these differences are that the Prandtl-Glauert ap- 
proximation, though valid for L^ach numbers in the neighborhood of 
unity, should be utilized only for very thin shapes; whereus the 
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von Karman method, though applicable to any reasonable "shape, is no 
longer valid for Mach numbers beyond the critical value. 

As final numerical applications of the results of this paper, 
the maximum values of the negative pressure coefficient *(^p,m)jj^3^* 

the critical pressure coefficient *^^p,m)q » limiting 

pressure coefficient ' (pl^^y^^n^ calculated for various 

values of the thickness coefficient and of the stream Mach number. 

The maximum values of • the negative pressure coefficient for 
various values of t and of M are obtained by means of equa* 
tions (17) and (20) • Table X shows values of a^^, and a^ 

defined, respectively, as the coefficients of t, t^, and t^ in 
equation (1?)- The corresponding values of Cp^|yf, obtained by 

means of equation (20), are listed in table XI, together v/ith values 

Cp^M calculated by the von Karman method. Figure I4. shows the 
variation of /C„ J\ with. Mach number for several values of t. 

The critical pressure coefficient (pp,M) calculated by 

means of the follov/ing expression obtained by substituting for q^ 
in equation (I9) the expression for q^^.^ from equation (I6): 



1 + 



2 + (y-1)m2 



(21) 



Table XII gives the values of ^^p,m} calculated by means of this 

equation, and figure h shov/s the corresponding granh. The inter- 
sections of the (^p,}J[J curves with the (^v^Ia) curve yield 

the critical values of the stream Mach number which are listed in 
table XII. 

As noted once before, for a given value of t, the series for 
q (equation (15)) apparently diverges for a definite value of the 
stream Mach number. It is reasonable to assume that this value of. 
the Mach number marks the limit of irrotational potential flov/ and 
also probably indicates, the first appearance of a compression shock 
at the solid boundary. Equation (I7) for <ij^Q^^ can be used to 

estimate the limiting values of M according to the follovvdng rough 
criterion. By means of table X, expressions for ^^x* 



form of power series- in t, can be obtained' f or any ^value of the 
stream Mach number in the range 0^ M <1." • Far" a given value of 
the thickness coefficient t, a series for q^^^^ can then be found* 
so that a term-f or-term comparison with the' harnionid series 

00 

V" 1 

y — yields a depreasing sequence for the ratio of corresponding 
n«l 

terms. The value of the Mach nuiriber thus -obtained* is chosen as the 
limiting value of the Mach number. TalDle XII lists both the values 

of (Cp m)t "^^^ corresponding values of M for different 

>■ ^ ' lim 

values of t. Figure I4. shows the curve connecting the limiting 

values of Cp^j^/^. The region between tVds curve and the (Cp,M) 

curve represents the supersonic range without compression shocks. 
It is emphasized that the mathematical procedure outlined in. this 
paragraph is a highly speculative one. 

The maxim'um speed that can be attained, by an adiabatic fluid 
is equal to the speed when p » p « c « 0 and is given by Bernoulli's 
equation 

2 = 1 



%hs = ^ * 



Values of the pressure coefficient 

are listed in table XIII for various values of the-'undisturbed- 

stream Mach number M, and figure I4. shows the corresponding curve. 
The region between this curve and the /c^ tvA curve reoresents 

the supersonic range with compression shock. On and beyond the 
(^P^m)^^^^ adiabatic fluid ceases to exist; that is^, ab- 

solute vacuum prevails. 

In conclusion, it may be remarked that the results of von Karma'n, 
shown by the small circles in figure I4., are obtained independently 
of any assumption concerning the shape of the solid boundary, whereas 
the results of this paper were • obtained for a specific family of 
shapes. As shown by the curves of figure k, nevertheless, the re- • 
suits of this paper agree with those obtained by von Karman/s method. 
This agreement has some justification, for the values of ^Cp 
and Mqj, depend mainly on the dimensions of a body - that is! 
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thickness coefficient - and not on its shape. The hodograph method 
as employed by von Karman, however, yields results that cease to be 
valid when the local Mach number equals or exceeds unity, whereas 
the iteration method utilized in the present paper is valid for local 
Mach numbers greater than unity and, for the family of shapes chosen, 
yields some information with regard to a supersonic region of flow 
without shock. The limiting value of M for a given shape, beyond 
v/hich supersonic flow without shock does not exist, appears to de- 
pend on the convergence of the power series in t for the velocity. 
Although only a few terms of the series have been obtained and 
therefore the limiting value of M cannot be given precisely, 
nevertheless it ,is believed that a reasonable estimate of the value 
for Miiui ^® made by the comparison test with the harmonic 

series. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va, 
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APPENDIX A 



THE INCOJiPRESSIBLE PLOViT PAST A CURVED SURFACE 



In the search for a shape v/hich satisfies the conditions that 
it be thin and that it possess no stagnation points, the first 
thought is of a straight-line segment. It is well knovm that a 
straight-line segment of length 40 is obtained from a circle of 
radius c by means of a Joiikowski transformation. If Z denotes 
the plane of the segment and Z' • the plane of the circle, then 

c2 

Z = Z» + — (A-1) 

The singular points of this transformation are determined by the 
equation 

dZ' \ z'/V z'/ 

which shows zeros at Z* » ±c. In order to raise the top surface 
and lower the bottom surface of the line segment, it is necessary 
only to place two additional zeros at Z* ±id where d c. 
Analogous to equation (A-^2), 



dZ 



H.(-F)f-Fr)(-^)(-F) 



Then, on integration of equation (A-3), it follows that 



o2 - d2 c2d2 



Z' + ^ + r (A-i+) 



The parametric equations of the shape in the Z plane corresponding to 
the circle of radius c in the Z' plane are obtained by substi- 
tuting Z' = oe-"-^ in equation (A-k) ; thus. 
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d2 

X = 2c cos 9 (3 cos 9 - cos 59) 

5o 



Y = — (3 sin e - sin 39) 

3c 



(A-5) 



The fsunily of shapes given by equations (A-5) includes, on one hand, 
the straight-line segment with d =» 0 and, on the other hand, a 
shape having four cusps syinnetrioally placed with respect to the co- 
ordinate axes with d » c. For 0< d <c, the slope dY/dX is 
zero for 6 =" 0, Tr/2, and ir. The shaoe thus has ousos at 

/ d2\ 

X « ±2c(l — — ], Y 0, and the ma-xiiaum and minimum points aire 



at X » 0, Y — , and 

3 c 



0, Y =» - — , resoec.tively. . 
3 c 



The complex potential for a circular cylinder of radius c, 
fixed in a stream of uniform velocity U in, the positive direc- 
tion of the real axis, is given by 




(A-6) 



The complex velocity past the corresponding shape in the Z plane is 



U - IV 



dZ 



dF dZ' 
dZ' dZ 



or 



u2 + » — — — — 
dZ» dZ' dZ dZ" 



By means of equations (A-i+) and (A-6), it follows that 
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q' 



,2 = u2 




and at the surface of the profile, where » ce^^, 

^ ■ — 72 — ^~7Ff 

1 + 2 — cos 29 '+ I — ] 

It v/ill be convenient to consider F, Z, and Z' as nondi- 

mensional quantities. Thus, in the plane Z' the unit of length 

is the radius c of the circle; whereas, in the plane Z, the unit 

of length is the semichord 's of the shape. Then Z*, Z, and P 

denote Z^/c, z/s, and f/Us, respectively. Also/ the ratio 

2/2 ' . 

d /c is designated by 

According to equations (A-5)# "the semichord s is given by 
and the thickness coefficient t by 

3 - e 



With the introduction of these nev^ designations, equations (A-ij.), 
(A-5), (A-6), and (A^?) become, respectively. 
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— : cos e - r (3 cos e - COS 36)1 

h h I 

} (A-9) 



Y = I (3 sin e - sin 39) 



J 



P=^^{Z'+— ) (A-10) 



2 1 /X 
« (A-11) 

1 + 2€ cos 2B + € 



As a nxmerical example of the use of equations (A-9)i (A-IO), 
and (A-11), table I gives the coordinates for the shape t » 0.10 
together with the velocity and pressure distributions along the pro- 
file/ The pressure coefficient Cp is calculated by means of 
Bernoulli's eqiaation: 

P - Po p 
|-PU2 

where 

p local static pressure 

p^ static pressure in undisturbed stream 

p density of fluid 

Figure 1 shov/s the curves of the shape aiid the velocity distribu- 
tion, and figure 5 shows the graph of the pressure distribution. 
Because the body is placed symmetrically with respect to the undis- 
turbed stream, the flow is identical v/ith that over a solid boundary 
composed of the X axis from X ^ to X = 1, the upper surface of 
the shape with 0 ^ 9 ^ tt, and the X axis from X =« -1 to X -a>. 
This boundary is called a bump. 
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APPENDIX B 



THE COMPRESSIBLE FLOW PAST A BUII> 



Before proceeding v/ith the iteration process, developments for 
and Y in positive, integral powers of the thickness parameter t 
must be obtained. From equation (A-8), 




(B-1) 



and from equation, (A-IO), 




(fl-2) 



By means of a Taylor expansion in the neighborhood of t « 0, it is 
possible to express P as a power series in t in which the coef- 
ficients are functions of Z. Thus, according to equation (B-l), 
t « 0 v/hen 




or 

Z» « Z + v^Z^ - i • 



where the positive sign is taken with the radical because the points 
at infinity of the Z and planes must correspond. Now ' 
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on the right-hand sides of equations (B-l) and (B-2), Z is con- 
sidered constant and Z » • is considered the oarameter, it follows 
that 



F = Z + 



iz.(z2-l)(z-yz2.l) 



t + . . . 



Since P = + iij; , virhere ^ is the streaia ftmction. 



^-x. it - (z2 - i)(z .#7"i )- (f - i)(z -vt^-i j 



t + . . . (B-3) 



In a similar manner, equations (A-9) can be written as follows: 
» 

X T COS 6 



t «• - 



sin^ 9 cos 9 



Y ? r (X - cos 9) tan 9 . 



If X is considered constant and 6 is considered the parameter, 
a Taylor expansion yields ' . 



y-(l-x2//^rt.3x2t2+^(8X^-3x2)t3 - . . 



(B-5) 



Determination of f^i 

By means of the transformation x = X, y ^ pY, where 
\A " M equation (7) for ff^ becomes 

6x2 . ^^2 . 



(B-6). 



By analogy with the coefficient of t in equation (3-3) for jzf, 
it is suggested. that 
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J2(l » k - (z2 , l)G - ^y^)- - - O] (B-7) 

where the coefficient k is determined by means of the boundary- 
condition 

6^ dY ^ 6^ , _ 

6x dX 6r 



or 



6x dx 6y 



(B-8) 



and where the boundary, obtained from equation (B>5), is given by 



y p(l - x2)5/^ [t - 3x2t2 . I (8x^ - Jx^) t?' . . 



(8-9) 



It is clear that the boundary condition need be satisfied only to 
the same degree in t as is involved in the development of the ve- 
locity potential ^. Thus, 



j^- :c + k jx - (z2 - l) (z - ^z2 -"i) - - i)(t - 7? - l) 



(z-n/z^-I )- 27(7-v^.l ) 



t + . . . (B-10) 



dfi ^ , i . 

Oz 08 '- 



+ N/a2l 1 G - v/z2 . 1 )+ 712-1 (f. i) 



t + ... (B-11) 



dy Adz 6F/ 

, \/i27I .. (z - yziTi )+ v^FTi (l , ^/?~l)|t + . . 



4^ » -3{Sx Vl - x2 t + . . . 
dx • 



. (B-12) 
(.B-13X. 
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Suppose that x = cos a. Then, insofar as the terms of the first 
power in t are concerned, it follows th&.% on the boundary, 

z » cos a + ipt sin^a 
F = cos a - ipt sin^a 

Vz^ - 1 = i sin a ^1 - ipt sin^a^ 

\/r^ - 1 =» -i sin a (l ipt sin^a) 
Then, on the boundary. 



r — « 1 - 5kt sin 2a 
ox 



dx 

^ « -5Ict sin 2a 
dy 

£l « -ifit sin 2a ■ 
dx 2 

The boimdary condition, equation (3-8), then yields 



2p 



Therefore 

2p L J 

This expression for 0 can be put into a simple form by means of 
the transformation 

z » cosh ^ (B-15) 

where 



2k 



Also, 

X =* cosh ^ cos T| ^ 

y » sinh sin T) J 
Then from equation (B-lU) it follows that 

JZ^Q = cosh cos T) 



* ^ cos T) - e"^^cos 3^1 



) 



(B-16) 



> (3-17) 



Inversion of Equations- (B-I6) 

The relationship between the rectangular Cartesian coordinates 
X, !l and the elliptic coordinates tj. is obtained as follows 

Invert equations (B-I6) and solve for and r|j tl\us 



^cosh §y 



and 



Vposh 



y 



vsinh 



Solve for sirJi , . 



2 sinh^£ =» -b + a/o^ + Uy^ 



and solve for sin^r) 



2 sin^n « b + Vb^ + l+y2 



where 



b « 1 - (x2 + y2) 
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By means of the transformation 

X = X 
y - ?Y 



J!j it follov/s that 



2 sinh^i = -b + + i;p2y2 
2 sin^rj = b + \^ + Up^Y^" I 



(B-18) 



v,rhere 

b = 1 - (x^ + p^Y^) 



26 



APPENDIX C 
DETERI-lIiJATION OP 



The differential equation (8) for j^fg terms of the co- 
ordinates X and y becomes 



6y2 



1 - f 



dx dx^ 



(C-l) 



By means of the symbolic relations 
6x hz hz 



+ 2 
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and the transformation 

z « cosh I 

equation (C-l) can be expressed in terms of the complex variables I 
and I as follov/s': 

sinh t sinh i^khi P 1 ^"""^ ^ 



+ 

sinh 



1 1 d^i^^i cosh t ejZ^i ^ 1 cosh ^ d)2^A 

^ ^F/ Vinh^t sinh^^ 6^. sinh^f cJp sinh? dl ^ 



\sinh ^ sinh [ 6l/\sinh2^ 6^^ 

cosh i hj^ 1 d% ^ cosht^6^j^ Y 

" sinh^t 6p ' sinh^f bV^ * sinh^t 6f /[ 

Now, equation (B-I7) for j^i can be rewritten as follows: 

The substitution of this expression for into the right-hand 

side of equation (C-2) then yields the following differential equa- 
tion for jzfg: 

ik.^.ll:^ . (3-,),^ e-3i^ - cos . 

+ (y-1)(i- P^)(e"5^ cos 3r\ - e'^cos (C-U) 
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The right-hand side of this equation suggests a solution of the form 
^ t) + Pj(£)cos 371 + p5(g)cos 5n (C-5) 

If this expression for inserted into equation (C-i4.) and the 

coefficients of cos "n, cos Jt), and cos 5*0 are equated on both 
sides of the equation, the following differential equations for 
Fj^(§), F^iQ, and FcjCS) obtained: 

^ 2 

, » iijli-. [(y+D ^. (5.Y)/|(e"^ - e"5|) (C-6) 
d§ • 8 ^ 

§ - 9F5 - I (r*l) (j-^j e-S , (0-7) 



2 2 
^ - 25F5 - - I (^)-"^ 

The solutions for F-^, Fz, and F^ are easily obtained and are 
as follows: 

Pl " ~ Ih^'- * Ry+D + (3-'Y)p2l (e-5S - 3e-^^)[- (c-9) 



2 n 



^2 
2 



^ 128 p4 

where A^, Aj, and A^ are arbitrary constants. It is noted 
that, in general, the expressions for P]_, Pj, . and F5 should 
each contain tvro arbitrary constants; however, the condition that 
Fi, F5, and F5 vanish at infinity requires that the omitted 

arbitrary constants be taken equal to zero. 
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The arbitrary constants A]^, Aj, and A^ are determined by 
means of the boundary equation (B-8), In terms of the variables ^ 
and 't] equation (3-8) takes the form 



sinh ^ cos r\ ^ - cosh § sin tq ^^5Z 
dE, 6r) / dx 



p^|cosh I sin ^1 §1 sinh g cos (C-12) 
and the velocity potential jl( has the form 
0 = cosh § cos r\ + — ^3©"*^ cos e cos 3^ 

Up 

+ t ^Fj^cos Tj + F^CQs 5^1 cos 5'n^+ • ..^ (C-13) 

With X =5 cos a, the boundary is given by 

y » pt sin5a - 3pt^sin\ cos^a + ^|3t^ ^Bsin^a cos^^a - 5sin5a cos^a^ - ♦ 
and 

£l =x ^Ipt sin 2a + ^?t^ (2sin 2a + 5sin Ua) 
dx 2^^ 8 



- -^pt^ (9sin Ua + lUsin 6a) + . • • (C-lij.) 
Id 



At the boundary, iiiclusive of terms containing the third power of t. 
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p = 1 - (x^ + y2) = sin^a (l - p^t^sin^a+ 6p^t^sin^a cos^a) 
^ = pt sin^a(l - 3t cos^a} 



Sin 



cosh ^ » 1 + ~ p^t^sin^a - 3p^t^sin\ cos^a 



r 2 1 P 2 ii 22 2 

a 1 _ pt sin a + p t sin a, + 3pt sin a cos a 



;> (C-15) 



sin "n » sm 



cos "0 =» cos 



a ^1 - I p^t^sin^a + i- p^t^sin^a 

+ jp^t^sin^a cos^a - 3p^t^sin^a oos^a^ 
° (■"■ " 2 P^^^s^^^'' * 5p^t^sin+a cos^o^ 



At the boundary then. 



dPi 3 1 - p2 I- 



2! 



T — s-Ai + U (y+1) + (3-y)P !f 



6U p' 



d£ 128 p^ L ^ ^ J 



dPq 9 1 - 



d^ 128 p^ 



[-5A5 - 3(y+i)(i-P^)] 



(c-16) 



Hence, with all terms involving powers of t higher than the second 
excluded, the boundary condition (equation (C-12)) yields the 
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following expressions for the arbitrary constants A]^,- Aj, and 



Ai = h FCy+D + (3-y)P^1 - 8 

t- J 1 - 



-^3 3.1,8^ 



^3 ^15 1 - p 



(c-17) 



The expressions for Pj^, F3, and then become: 



Pl(g) - - |(y+1) + (3-Y)P^l(iia'^ - 30"^^ + 

614. - 



^ 128 



{Y^l)(i^j(3o-5^-5e-59 



. 1. LUJllg^ e-5S 
16 p2 



Fail) " -2- (y+1) 

6ij.O 



. 5 3 3P ^- -5S 
+ — A e 

80 



> (C-18) 
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Equation (C-13), together with equations (C-18), represents a solu- 
tion of the fundamental differential equation (1) and satisfies the 
boundary conditions at the surface of the solid boundary and at 
infinity, insofar as the terms inclusive of the second pov/er in t 
are concerned. The present process can be extended to the higher 
terms in the development of the velocity potential 0, but it can 
readily be seen from the complexity of the right-hand side of equa- 
tion (9) that the labor involved would increase rapidly with the 
degree of approximation. In viev/ of the importance - of the problem 
it has nevertheless been thought worth while to extend the calcu- 
lation to the third approximation ^-z* 
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APPENDIX D 

DSTERMINATIOIf- OF- 0^ ' 

The differential equation (9) for 0^ in terms, of the complex 
variables z and z can be written as follows: 

* (i^iz - )^iz)02^2zz - i2.r\ \ ......... ^^"^^ 

* >• ... * 

where, for example. 



^52 2* 



3 



Again, complex variables and t are introduced in place 
of z and z by means of the transformation 



z « c 



osh % 



3h 



According to equation (C-3), 



(Cr.3) 



Consider the term 



Now 



5 e-2t 



■ lZ2 



2p sinh i 



and 



8f 



^Izi^lz " ^'^^ 
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^ _ ^ e^^sin 3*0 - e"^^ sin t) 

2p sinh ^ sinh J 



It follows then that 



Ti sinh ^ sinh J 



+ ^(Y+l)(l-p^)e"5S -1(y+1)(i-P^) + SP^J e*?!! cos 3t1 
- (Y+l)(l-p2)^e"7S cos 5-0 +(y+1)(i-P^)^®~^^ cos Tn 



(D-2) 



In the handling of the terms on the right-hand side of equa- 
tion (D-l) that contain the derivatives of ^^.f is convenient 
to separate ^2 into two parts: namely, ^2.^ » which is a function 
of t_plus a function of t, sind which is a function of ^ 

and L Thus, if the variables S and t are introduced into 
the expression for follows easily that 

JzTp" -c(l^-^3e-^-2r^^.2t-3l,i^-r.3e-2M,3-3t^l) L (^,.3) 
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32 



l6p 



9 
1280 



\2 



160^*= 



128 



\ f / 32 r 



D « -2- (y+1) (^-^ 

256 \ p'^ 



Table II gives values of A, B, C, and D for values of 
the Mach number M ranging from zero to nearly unity. The value 
of Y takan is l.i4.05« 

Consider the term 
T2 - [(y+1)(w2) 2p2] [(/u,*^i2^)(p<'2'. + i4'z) ' 

- \{f^lzz ' /lzz)(/2's- ^z'zY {hz -^l^X^2'zz - ^2'zl)\ 
+ 2 j( Y+1 ) (l+p2) - 2p2 J (^^^ ^ pTj^^ ) 
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Now 



2 zz 



sinh t 



and 



•* " ISiAe"^^ sin 2t) + 60iBe*^^r,in Urj 



sinh I sinh t 



+ 120B ^e'^^^cos 5ti - e'^^cos Jti^J 



=rl2A(e^^sin 3^] - e'^^sin r]) 



sinh ^ sinti ^ 



1208(0*"^^ sin 5ti - 5^^^ 



It follows that 
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T2 sinh t sinh t - ^ y [(y+1)(i-?^) + hf; j Ae'^S 

■ _ ^ - "j 



H. <|^[^Y+l)(l-^2)+Up^]B-||(Y+l)(l+P^)-2p2j(3C+2D)^^^^ 



||(r+l)(l+P^) - (5C + 2D) > 



+• J- iI5 (Y^i)(i_p2\3^ _2p2jD|e-^ cos 7^1 (D-J+) 



cos 5'n 



The calculation of *fx '^2 heen very involved. 

The calculations for the remaining terms on the right-hand side of 
equation (D-l), however, are quite laborious and therefore only the 
final results v/ill be presented. The terms to be oalcvilated are 
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The expressions for the constants AqP are as follov/s: 
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|(y+1)(i-P^)^ + I + -(Y+l)(l-p2)(6A-15B+28C -7D) 



+ 12p (3y + 22)C + 6a - 15B - 8d] 
A 1 « . (y+1)(i-P^) + I (y+1)(i-P^)(30B - 19C - 2i|D) 



- 12p |^y(3C + 2D) 30B + 7C + 20dJ 



i5 



(y+1)(i-P^)(3C - 2D) - 120p(2C - D) 



Ai3^ =^ [(Y+i)(i-p2) + up2]c ; 

iJ^ -^ (Y+l)(l-pi+) + 2 (.^l)(l-p2)(2A - 103 + 5C) 

^ 32p^ P 

+ 12p j_Y(3C + 2D) - 30B - 5C + 12d| • 

V ^ i(Y+l)(l-p2f ^ 8p2]+ 1 (Y+l)(l-p2)(l5B - 16C - D) 

6i+p " P 

+ 12p ['(-3y + 10)C + 15bJ 
^11^ ' - y (y+1)(i-P^)(3C - 2D) - l68p(2C - D) 

Ai5^ - |^[(y+1)(i-P^) + i^P^lc 



50- 



V ' - I (Y+l)(l-i3^)(5A + iiC - 6D) 

a/ = - ^ 1^ (Y+l)(l-p2)(9B - D) 

- 12pii2C + (2y - 23)d] 
A^5 = _ ^ (Y+l)(l-p2)(C - 2D) + 252pC 

T [(y*i)(i-P^) + hp^]ic - 2D) 
^1/ * - ^ (r+l)(l-p^)(3C - 2D) - 2l6p(2C - D) 



.7 - 27 



A5' « (Y+l)(l-p2) ^ 2 (Y+l)(l-p2V[^5B + 12C-29D) + 12p(2Y + l)D 

Ag^ - - .192p(C - 2D) 



A 7 60 

^13 " y 



_(Y+l)(l-p2) + i/](C - 2D) 



66 r 
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APPEiroiX E 



THE FLOW PAST A CUR\TED SURFACE BY pE POGGI I.ETHOD 



Equation (D-21) represents a solution of the fundamental dif- 
ferentie.1 equation (l) that satisfies the boundary conditions at the 
surface of the solid body and at infinity, insofar as the terms 
inclusive of the third pov/er of the . thickness coefficient t are 
concerned. The method used has been called Ackeret's iteration 
process and is valid for all values of the I^ch number M from zero 
to ^anity. On the other hand, the method of Poggi yields the com- 
ponents of the fluid velocity in the form of povrer series in I^. 
Since both Ackeret's and Poggi *s methods provide solutions of equa- 
tion (l), the two solutions must agree in the region common to both. 
The flow past the shape treated in the present paper will be calcu- 
lated by means of Poggi 's method and compared v/ith that obtained by 
Ackeret^s method. 



Poggi 's method consists in regarding a compressible fluid as 
an incompressible fluid with a continuous distribution of sources 
in the region external to the solid boundary. In orvder to express 
the intensity of the source distribution, it is first necessary to 
determine the incompressible flow, vv'hich serves as the zero approxi- 
mation to the flow of the compressible fluid. The first-order 
effect of compressibility on the velocity of the fluid is then given 
by a set of double integrals extended over the entire region of flow. 
In reference 5 "^he surface integrals are replaced by line integrals, 
which are evaluated by the methods of the calculus of residues. 
For the example treated herein, the general results given in refer- 
ence 5 are inunediately applicable and are as follov/s: 

Let Z - Z(Z*) be the conformal transformation of the profile 
in the Z plane into the circle of radius c with center at the 
origin of the Z' plane. Then, if V/^ is the complex velocity of 
an incompressible fluid past the circle, the following expressions 
are form.ed: 



dZ dZ dZ' \ dZ / 
dZ dZ dZ' \ dz/ J 



!!4L 

° dZ 



> 



Wo^^dZ. 



(E-1) 



52 



The complex velocity V( of the compressible fluid past the profile 
in the Z plane is given by 



(E-2) 



where only the terms involving the square of the Mach number 
shovm. The expression for Wj^ is as follovz-s: 



are 



W]_ « — 



2 f 2 



(E-3) 



where F| — ,Z'I is obbaiued from F(Z',20 by replacing Z» by 

c2 1 / c2\ 
~; S(Z') denotes the sum of tlie residues of = P(z* ) 

at the poles v^ithin the circle C]_ of radius c, the subscript P 



being dropped after the evaluation; and S 



is obtained from 



the expression for S(Z') by replacing Z* by . 

2' 

For the present example, according to equations (A-3)# (A-U), 
and (A-6), 



dz' V 2'y\ z'./\ z'A zv 



o2. - d2 
2 » Z» + • + 

Z» 



c a 
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,5 



(A-3) 

(A-U) 



and 



51; 

Considsr noxv the contour integral 



1 f(z»,— )dZ« = 2TTiS(2»p) 



fCi 

In the first term on the right-hand side of equation (E-5) only the 
simple poles et Z' =» ±icl are internal to the circle C-j^* 
contributions of these simple poles to the residue are, at Z* = id, 

U^c^ 1 



2id(c2 - d^) Z-p - id 



and, at Z* =» -?id, 

u3c^ 



2id(c2 - d^) Z»p + id 



If the two expressions are added and the subscript P is dropped, 
the result is 

^''"^ ' (E-6) 

The second tarm on the right-hand side of equation (E-5) can be re-* 
vo'itten as follows: 

2«' # --iu5d(o2* f XOS--^ 

(Z' - id)'^(Z' + idr (Z'- id)'^(Z' + id)'' il-+z 

d • 



The contributions of the double poles at Z' = id and Z' = -id 
to the residue are, respectively. 
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and 



U^c^ Z'p + 2id 1 iU^c^d i 



2id (Z«^ + idr 2 c2 _ Z'p + id 



c2 - d2 



^ .(Z.p + id)2 %2^^2 

If the two expressions are added and the subscript P is dropped, 
the result is 

au^c^d^ IT^cSd^ 1 • 



(s''- + d'^) C- - d" Z'^- + d^ 



;.2-d2 c2-d2 



+ d^; + d^ 

If the expressions (B-6) and (E-7) are added, it folloxvs that 

S(Z')-Uc fx* i-) / ^ on2 (E-8) 

From equation (E-i).), 

z.2(z.2.d^)(z.a*d2) 



If S(*Z») is formed from equation (E-8) and Z"* is replaced by 
— , it follows that 
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c2*d2, c2 -d^^^■o''-d22■2)z.2 

lo/r I- 



2c2 



(E-10) 



Finally, if P(Z',Z«) is fox'med fi'om equation (E-9) and Z' is 
replaced by e^/Z', it follws that 



F 



log 



Z» + id 



(c^ + d^Z'^f \ 2id, Z' - id 



(E-11) 



By means of equations (E-o), (E-9), (E-IO), and (B-ll) and 
placing Z«/c, Z'/c, V/j^/u, and d^/c^ by Z', Z', Wj_, 
respectively, it folloxvs that 

/' ■ . ' 

it IV 2 1- £7 (2,2. ,)2 (■„,2,^^^ 



by re- 
and c , 



7.2 . 



2e 



(c + Z'2)(ez»2 + l) (eZ'^ + l)' 



Z'2 + (1 + €} 



s Z.» 



2 ^ + 

2iVc " Z« - i\/e 



los 



. 1 + -i- + z^^(z«s - l) 

(Z'2 .c)(z.2 . c) . 



2€ 



Z.Z' ^ (1 ^c) log 



. 2iv^ Z» - iv€ 



(E-12) 



The complex velocity Yi of the compressible fluid past the 
profile in the Z plane is given by equation (E-2) v/here 



dZ« 



Z 



dZ . (z«^ - OC''''^ * 0 



(from equation (A-5)) (E-IJ) 
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Wo - 1 - 



(from equation (A-6)) (E-II4.) 



and Wj^ is given by equation (E-12). 

The velocity at the surface of the profile is obtained; by. means 
of equations (E-2), (E-12), (E-13), and (E-li+) v/ith Z' » e^ ■ and 
Z"' " e'^^ . Thus 



^profile 



1 . 

v( + 2cco3 29^?.V l + 2e COS 



X + 



€(1 - € - 2cos 29) 



(1 * 



1 - 



2(1-0: 



1 + 2e cos .29 + € 



log 



1 + € 
1 - € 



(i>0 • 1 

2>/i' l + 2e CO."; 29 + C*^ 



2e. + (i + c^) 

COS 29 - 1 + C +2V^ sin 9 
— log ' ^^ ' — 7 

sin 6 • -l + c rSV^ ' sin 9 



1/1 ^\ D -I- -1 2\/c cos 9 
- kd - r) coso tan -X.^ . — ^ 

1 - € 



. ■•■ CE-15) 



For 6 » 0, 



.\2 



or .1 



log 



1 + € ■ 
1 - € 



■ 2V€ [ 



(E-l6) 
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For e = 



- 1 + 



_ £1jli1 log ^ ^' 

(1 - €)2 1 > ^ 



The velocity, along the. X axis external to the profile is 
obtained from' equations (E-2), (E-12), (E-13), and (E-llj.) by taking 
Z;« = 2« = X'. Thus 



X' 



- 1 + 



x»' 



2c(l - €2)X>3(x»2 ^ 1) / i_j_c ,1 2vr X' 
- : / X • + tan 



2 

X'^ - € 




(E-18) 



The velocity along the Y axis external to the arofile is 
obtained .from equations (E-2.), (E-12), (E-13), and (E-UO by taking 
Z' » iy and- Z' « -lY'. Thus 
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Y' 



-.2. 




V 2 l-c/ (c - Y.2)'^(l , €Y.2y 



Y.2- 



- 1 + 



(c - Y'2)(l -cY'2) . (c , y.2y 



Y' -Vi:y (y2 + i)(i . eY'2)(c . y'2)' 



f 



(E-19) 



Table III gives values of the velocity corresponding to points along 

the profile for the nwnericai case e= i (or t = O.IO) at Jfeoh 

numbers M =0.50 and O.75. Figure 1 shows a comparison of the 
velocity distribution for M = O.50 with the corresponding calcu- 
lation according to the iteration method. 

It is noted that by Poggi's method the solution of the problem 
is given by the components of the fluid velocity, whereas bv' i^ckeret's 
method the solution is obtained in the form of the velocitv" ootential 
of the fluid motion. Before a comparison of the t^vo methods can be 
made therefore, it is necessary to obtain from the velocity potential 
^ol Ackeret; the velocity components along the coordinate axes. 
This calculation is performed in appendix F. 
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APPEiroiX F 
DETERMIKATIOH OF TSB VELOCITY COiffOlENTS 

In terms of the complex variables t and t, the velocity 
components in the direction of the coordinate axes are given by 

' ^ ._l ^ 1. ¥ 



u » 



V « — » ip 



oX sinh i bt sinh 1 6 1 

)/ 



(P-1) 



where the velocity potential jzf is obtained from equation (D-21), 
According to equations (6) 



u » I + till + t^Ug + t-^Uj + 

V = tV|^ + f^v^ + t-^v.j + . . 



I 



where, in general. 



6j< 



n 



.6x 



(F-2) 



In terms of the complex variables t and 

.i(3e-t..-3^V^.e-^) 



Therefore, 



5 -22. 
Ut - — e ^ cos 2ti 

2p 



3 -2^ 
V, o - e ^ sin 2ri 



> 



(F-3) 



Similarly, with reference to equations (p-3), it follows that: 
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+ 3Ae"2S +i4i)e"*'^)co3 2ti + 1+11406 "^^ + (D - 15B)e"^S]cos l+rj 

I** — 

52 sinh% cosh 2D * c(2e''^^ + e"^^) 
cosh 2g - cos 211 



V2 » i+p sin 2ii l(liD - 3Ae-^^ + ItDe'^''^) 



+ 2 



UDe"^^ - (D + 15B 



COS 2t) 



1+ sinh%i2D V c(2e-2? + e'^^)] 



cosh 2£ - cos 2r| 



■J 



The general expressions for the velocity components and 
v-z are too cumber some to be given here* Instead, only the expres- 
sions for the velocity components along the profile vail be given. 
Thus, along the profile/ if powers of t higher than the third are 
neglected, 

N ■ 

« (^i + 5G3 + 5G5 ^ 7G^)^ + 2(3G^ + 7Gy^^ cos .2a 



(F-5) 



V5 « i-2ti j(^Gj» + + Gj^) sin 2a 



where the primes denote differentiation with regard to the inde- 
pendent variable £, and the zero subscripts denote evaluation for 
g « 0, Explicit expressions for (^Gj^^q* (^^)o* (^5}o' (^7' ' 

(^lOo' (^3')o' (^50o' (^7*^0 ^* 

this appendix. 
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In order to compare the velocity of the fluid given by equa- 
tions (P-6) with that obtained by Poggi's method given by equa- 
tion (E-I5), explicit expressions for equations (F-6), v/ith powers 
of the Mach nuraber higher than the second neglected, are obtained 
and may be expressed as follows: 



cos 2a 



\ h 16 J \ 2 h 32/ 

H. (^t^ , ^ t^^; COS ha - 201 ^3 
V 16 / 32 

. t2 - ill tA . (If t . 5 ,2 . 369 ,3) ,,3 2a 

,('5it2 .63,3\,,3i^^.|j5,3 6a 
\16 32 / 6U 

V = t . ^ t2 . 1 t3)sin 2a . ^3t2 -g t3)sin Ua - ^ t3 sin 6a 



)(F- 



t3 sin 2a + ( -2- , 2 t3^ sin i;a 2 t3 sin 6a 

32 Vi6 ii y i+ 



The magnitude q of the velocity, to the same degree of approxi- 
mation, is then given by 

- ^ ^ (-1 H^'^'i - ^ (i -I ^ 

?5 t2 . m t3) . (1 i t . i t2 - 873 t5 Vos 2a 
Vs 160 / \ 8 6itO / 



- t3 cos 6a + 

6k 



I. — t ^ * ) COS i+a =^ t-' cos 6a 

\16 32 / 128 



(F-6) 



v/hich is to be compared with the corresponding expression obtained 
by Poggi's method. Thus, if in equation (E-I5) € is replaced by 
3tX2 + t) and all terms are expanded in powers of t up to and 
including t^, the resulting equa-tion is as follows; 



6]+ 



^profile = * l6 * - i6 V V2 * IT * IT y 



cos 
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16 ^ ^ 6U 
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^27Y2.t3) COSI46 -Ii2it5 cos 6 
16 ^ V 128 
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(P-9) 



From the first of equations (A-9) with X =» cos a, it follows that 

COS a ( 1 - t sin^Q) cos Q • 
and therefore . . 

cos 2a « cos 28 - sin229 - i- t2 sin2e^ 

cos Ua =» cos ii.9 - 2t sin UG sin 2B 
cos 6a =» cos 69 

If cos 2a, cos ha, and cos 6a in equation (F-8) are replaced by 
these expressions, equation (F-'9) again obtained. 

The expressions for (O^)^, (g^)^, (65)^, (g^)^, (G^t)^, 
(*^30 ' C^5')- ' ^^^^ (^7*^0 follov/s: 



0 " 32 



2^ P 2p^ 2 8 



(y+1) y(^-|^*Tf B-|c + |Dy+ ti(-6A + 15B-Ic + ll4D^ 
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11 5 1 13 



16 



2\2 

«2 



256p- 



(y+i)(i-p^)-^(y+i) 



b3 80p '*^uo 



+ (y+1) — (fA-fB+ — C 



+ P^-9B + i- C + 2D^ 



5 


£5 


160 


2 



P 2p^ p2 2 \p2^ 



'0 1 



1^ 



9 
1280 



(y+i) + (y+1) a + 2.b*-Id\ , i 

5 P \ 20 i; 20 / 5 



(^7)0 - i 



4p.^.i.J^.iL.|(,.i(4tp.^(Y+i/^^^ 

S V/ 32 P 



8 ^ P 8p5 I^p2 





9 


p2 


1792 



(y+1) 4+1 YPD-;^ (y+1) (U5B - 17D) * ^ PD 
p-' 7 28 p 7 



V 1 /o 32 



. 5 , . , . 1 . ^ . ^(Y+i)f'4lp * i(Y+/ 4t ^ 



P 2p5 2 



^P 



8 



(G3O =-2- 
^ ^ ^0 52 



11 



1 3 



~ P + 5 - ~ - - ^(r+i) -r P - ^(y+1) 



p2 2p3 2 Vpi 



8 



'it 1 
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TABLE I 



VELOCITY AND PRESSURE DISTRIBUTIONS AT THE SURFACE 
OF A BULIP, t = 0.10 



e 

(deg) 


X 

/•n j_ »■ t k r\\\ 

(Equation (A-9) ) 


1 

(Equation (A-9)) 


q 

(Equation (A-IO)) 


Cp 

(Equation (A-ll)) 


0 


1.000 


0 


.3750 


.2344 


5 


.9954 


.0001 


'.8764 


.2318 


10 


.9818 


.0005 


.8808 


.2241 


15 


.9595 


.0017 


.8881. 


.2113 


20 


.9287 


. .0040 


.8983 


.1931 


30 


.8444 


.0125 


.9272 


.1404 


40 


.7344 


-.0266 


.9667 


.0654 


50 


.6051 


.0450 


:1.014 • 


-.0292 


60 


.4625 


.0650 


■ '1.068 


-.1395 


70 


.3118 


.0830 


1.117 


-.2476 


80 


.1568 


.0955 


1.153 


-.3299 


90 


0 


.1000 

• 


1.167 

... 


-.3611 



TABLE II 

Vj-iLUES FOR A, B, C, AND D OBTAIIJED FROM EQUATION (D-3) 



M 


P 


A 


B 


C 


D 


0 


1 


0.06250 


-0.06250 


0 


0 


.10 


.99499 


.06221 


-.06285 


.00095 


0 


.20 


.97930 


.06160 


-.06390 


.00400 


.00015 


.30 


.95594 


.06169 


-.05570 


.00982 


.00083 


.40 


.91652 


.064S9 


-.06830 


.01990 


.00307 


.50 


.86603 


.07788 


-..07171 


.03751 


.00939 


.60 


.80000 


.11322 


-.07551 


.07057 


.02675 


.70 


.71414 


.24561 


-.07741 


.14211 


.07805 


.75 


.65144 


.40177 


-.07394 


.21371 


.13977 


.80 


.60000 


.73.343 


-.06322 


.34482 


.26722 


.83 


.55776 


1.11855 


-.04583 


.48400 


.41460 


.85 


.52678 


1.53475 


-.02476 


.62619 


.57315 


.90 


.43589 


4.07940 


.12446 


1.4241 


1.5367 


.92 


.39192 


6.81182 


.29917 


2.2282 


2.5674 


.94 


.34117 


12.9343 


,70926 


3.9598 


4.8722 


.96 


.28000 


31.0576 


1.97342 


8.8910 


11.683 


.93 


.19900 


152.428 


9.30328 


34.423 


49.732 



TABLE III 
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VELOCITY DISTRISUTIOIT AT TIE SURFACE OF A BUMP, t = 0.10, 
. .. ACCORDING TO T;iS POGGI IIETHOD 



9 

(flop;) 


i 




X 




,q 

(incom- 
pressible) 


Coeffi- 
cient 
of m2 


q 

(compressible) 
(Equation (E-17)) 


M = 0.50 


14 = 0.75 


0 


1.000 


0 


1.000 


0 


0.8750 


-0.0540 


0.861i5 


0.8446 


5 


,9962 


.0872 


.9954 


.0001 


.8765 


-.0535 


.8631 


.8464 


10 


.9848 


.1737 


.9818 


.0005 


.8808 


-.0521 


.8678 


.8515 


15 


.9659 


.2588 


.9595 


.0017 


.8881 


-.0497 


. 8757 


.8601 


20 


.9397 


.3420 


.9287 


.0040 


.3985 


-.0462 


.8867 


.8723 


30 


.8660 


.5000 


.8444 


.0125 


..9272 


-.0359 


.9182 


.9070 


40 


.7660 


.6423 


.7344 


.0266 


.9667 


-.0192. 


.9619 


.9559 


50 . 


.6428 


.7660 


.6051 


.0450 


•1.015 


.0039 


1.016 


1.017 


60 


.5000 


• .8660 


.4625 


.0650 


1.068 


.0330 


1.076 . 


1.086 


70 


.3420 


.9397 


.3118 


,0830 


1.117 


.0545 


1.133 


1.153 


80 


.1737 


.9848 


.1568 


.0955 


1.153 


.0795 


1.173 


1.198 


90 


0 


1.000 


0 


.1000 


1.167 


.1002 


1.192 


1.223 



■ ■ • TABLE ' lY ■ ■ ' ' . 

V/J.UES OF (Gi)q, (G3)o, (Gg)^, AND- (Gy)^ 
G Pmi AT END OF ^VITEITOIX P 



M 




«'1>0 


■ («3)o 


, /^.5)o 


(^7)o 


0.50 . 


■ 0.86603 


0.10344 


-0.79011 


. 1.6714 


-0.91449 


• .75 . 


.65144 


.2.3075 


-2.1109' 


2.8975 


-1.9512 


.83 


.55776 


.9.8911 • 


-6.7071 


5.9329 


-4.0659 


.90 . 


.43589 


68,849 


.-42,516 


26. '227 


-19.549 
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VALUES OF &i, ag," AJfP . OBTAIJIED FROM EQUATIONS (14) 





•• 0.50 


0.75 


0.83 


0.90 


ai 


0 


1.7321 


' 2.3026 


2.6893 


3.4412 


.1 


1.6974 


2.2565 


2.6355 


3.3724 


,2 


1.5935 


2.1184 : 


2.4742 


3.1659 


.3 , 


1.4203 


• 1.3881 


2.2052 


2.8218 


.4 


. 1.1778 


1.5658 


1.8287 


2.3400 


.5 


. .86603 


1.1513 


■ .1.3447- .. 


1.7206 


. 0 


.484S7 


. .64473 


.75501 


.96355 


.7 • 


.03464 


.04605 


.05379 


.06882 


.8 


•• -.48497- 


-.64473 


-.75301 


-.96355 


.9 


.-1.0739 • 


-1.4273 


-1.6674 


-2.1336 


.975 


-1.5610 


-2.0752 


■ -2.4237 


-3,1014 


1.0 


-1.7321 


-2.3026 


-2,6893 


-3.4412 


• as 


0 . 


2*2743 . 


6.4350 


10.136 


19.268 


.1 


2.0399 


4.9892 


9.3S21 


18.076 


.2 


1.3673 


3.7088 


7.2538 


14.643 


.3 


.34809 


1.7641 


4.0015 


9.4302 


.4 


-.86531 


-.56063 


.10228 


3.1662 


.5 


-2.0593 


-2.3678 


-3.7900 


-3.1035 


.6 


-2.9594 


-4.6459 


-6.8342 


-8.0405 


.7 


-3.2325 


-5.2746 


-8.0105 


-10.022 


.3 


-2.4744 


-4.0011 


-6.0811 


-7.0773 


.9 


-.23520 


.01270 


.33212 


2.9851 


.975 


2.7276 


5.4013 


9.0381 


16.723 


1.0 


' 4.0063 


7.7376 


. 12.825 


22.703 


as 


0 


3i4379 


21.877 


72.397 


4i57.84 


.1 


2.5124 


18.507 


63.3-33 


386.94 


.2 


,07542 


9,4933 


58.871 


251.93 


.3 


-2.9313 


-2.1171 


6.6216 


72.788 


..4 


-6.1915 


-12.039 


-22.646 


-92.533 


.5 


-5.3952 


-15.937 


-37.804 


-184.67 


.6 


-2.7710 


-11.097 


-31.791 


-164.59 


.7 


2.2311 


1.4456 


-5.9844 


-41.555 


.8 


7.0915 


14.256 


22.646 


92.580 


.9 


6.2552 


10.459 


14.506 


23.206 


.975 


-.82172 


-20.694 


-56.978 


-385.70 


1.0 


-10.033 


-39,831 


-101.38 


' -633.92 



TABLE VI 



VELOCITY DISTRIBUTION FOR A BUIiP., t = 0.10, CALCULATED 
BY imkllS OF TABLE V AITO EQUATION (13) 



q 





0.50 


0.75 


0.83 


0 


1.199 


1.307 


1.443 


.1 


1.193 


1.294 


1,421 


.2 


1.173 


1.258. 


1.359 


.3 


1.14.3 


1.204 


1.267 


.4 - 


1.104 


1.139 


1.159 


.5 


1.061 


1.071 


1.059 


.6 


1.016 


1.007 


.9752 


.7 


.9734 


.9533 


.9193 


.,8 


.9339 


.9098 


.8865 


.9 


.8965 


.8676 


. 8444 


.975 


.8704 


.8258 


.7910 


1.0 


. 8568 


.8072 


.7579 



• TABLE VII 

VALUES OF CRITICAL VELOCITY OBTAIIIED FROM EQUATION (16) 
j\m MjIXIMM VELOCITIES FOR A BULIP, t = 0.10, 
OBTAINED PROM EQUATIONS (17) .1ND (E-17) 



M 



0 

.2 
.3 
.4 
.5 
• 6 
.7 
.8 
.85 
.90 
,95 
1.0 



Iteration 
nothoc' 



1.137 
1.171 
1.176 
1.185 
1.199 
1.222 
1.265 
1.371 
1.523 



Poggi 
not hod 



1.167 
1.171 
1.176 
1.183 
1.192 
1.203 
1.216 
1.231 
1.239 



icr 



4.578 

3.067 
.2.516 
1.869 
1.574 
1.366 
1.212 
1.145 
1.093 
1.044 
1.000 
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VALUES OF q/qcr 



TABLR VIII 
FOR A BUWIP, t 



= 0.10, VaTH U = 0.83 



X 


q 


q/qcr 


0 


1.443 


1.230 


.1 


1.421 


1.211 




1 .359 


JL • low 


.3 


1.267 


1.080 


.4 


1.159 


f9881 


.5 


1.059 


.9028 


.6 


.9752 


.8314. 


.7 


.9193 


.7837 


.8 


.8665 


.7558 


.9 


. 8444 


*7200 


.975 


,7910 


.6744 


1.0 


.7579 


.646\ 



T/vBLE IX 



COM'.'lRISOK OF TIE PRESSURE DISTRIBUTION AT TIE SURFACE OP A BUMP, t = 0.10 
FOR M = 0,83 0BTAIHI3D BY MEANS , OP ^ TIE ITERATION, THE PRANDTL-GLAUBRT , 

Airo TIE VC/J liAPJ/IAN METHODS 



X 


Cp,H 


Cp,0 


Iteration 
method 
(Equation (20)) 


Prandtl- 
Glauert 
method 


von I£annan 
method 


0 


-0.9133 


-0,6342 


-0,7376 


-0.3537 


.1 


-.8677 


-.6142 


-.7107 


-.3426 


.2 


-.7389 


-.5558 


-.6337 


-.3100 


.5 


-.5480 


-.4665 , 


-.5202 


-.2602 


.4 


-.3294 


-.3514 


-.3810 


-.1960 


.5 


-.1134 


-.2150 


-.2257 


-.1199 


.6 


.0514 


-.0647 


-.0656 


-.0361 


.7 


.1617 


.0765 


.0752 


.0426 


.8 


.2221 


.2062 • 


.1972 


,1150 


.9 


.2868 


• ,3196 , 


.2983 . 


.1783 


.975 


.4043 


.3935; 


.3620 


.2195 


1.0 


,4708 • 


.4163 


. .3812 - 


.2322 



TABLE X 

VALUES OF a^, ag,.. Aim ag CALCULATED FROM EQUATION (17) 



. M . ■ 


. «1 ■ 


a2 


as 


0 


1.50000 


1.50000 


1 . 50000 


.2 


1.53093 


1,58726 


1.67024 


.3 


1.57243 


1.71149 


1,95357 


.4- 


1.63663 


1,92091 


2.43349 


.5 


1.73205 


. 2.27425 


3.49200 


.6 


1.87500 


2.91036 


5 . 75400. 


.7 


2 i 10042 


4.22402 


12.56611 


.8 


2 . 50000 


• 7.76704 


43.7829 


.85 


2.84747 ■ 


12.52960 


112.7760 . 


.90 


3.44124 


. 25.7418 


446.847 


.92 


3.82733 


. 38.3963 


964.722 



TiVBLE XI 



tlAXIMUlJl" VALUES OF THE PRESSUiRE COEFFICIENT Cp^M CALCULATED 
BY ISEMIS OF EQUATION (20) 



-(Cn m) 


\ 


0.05 


0,05 


0.08 


0.10 


.0,12 


0,12 


0.15 






(a) 








(a) 




0 


0.16406 


Oi 16406 


0.27744 


0.36000 


0,44856 


0,44856 


0,59344 


.2 


.16764 


.16773 


. 28369 


.36830 


.45913 


.45994 


.60787 


.3 


.17243 


.17266 


.29221 


.37964 


,47363 


.47537 


.62779 


.4 .. 


.18008 


,18035 


.30570 


.39771 


,49688 


.49933 


.66002 


.5 


.19166 


.19188 


.32661 


,42506 


.53375 


,53657 


,71196 


.6 


.20971 


. .20937 


.36000 


.47200 


..59437 


.59401 


.79907 


.7 


. .24055 


.23753 


.41994 


,55715 


.71014 


.69006 


.97236 


.8 


.30606 


.28925 


.56154 


.77056 


1.0165 


.87903 


1,4634 


.85 


.38190 


.35321 


,74868 


1.0737 


•1,4759 


1,0664 


2,2485 


.90 


..59091 


.42108 


1,3551 


2.1252 










• .85480 


.47966 












\ t 


0.18 


0.18 


0.20 


0.21 


0.22 


0,25. 


0.25 






(a) 








.(a) 


0 


0.75384 


0.75384 


•0.07000 


0.93098 


0.99396 


1.1953; 


1.1953 


O 


. .77275 


.77541 


' .89226 


,95504 


1,0199 


1,2274 


1.2352 


.3 


.79399 


■ .80489 


.92326 


.98859 


1,0561 


1.2724 . 


1.2903 


. .4 


.84179 


• .85192 


.97408 


1.0437 


1.1158' 


1 . 3470 


1,3793 


.5 


.91186 


.92436 


1,0581 


1,1353 


1,2154 


1.4732 


1.5208 


• .6 


1,0315 


1.0403 


1,2031 


1,2943 


1,3891 


1.6964 


1,7566 


.7 


1.2782 


1.2431 


1.5088 


1,6326 


1,7623 


2.1889 


2.1973 


.8 • 


2.0175 


1.6781 


2.4539 


2,6938 


2,9491 


3.8118 

1 


3,3117 



^Method of von Ka^rnufn. 
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TABLE XII 



CRITICAL Airo LIMITING VALUES OF M /J-ID CRITICAL V/iLUES OP Cp^M 
CALCULATED BY MEANS OF -EQUATION (21) 



M 


-(Cp,M)cr 


t 


i. 9 


Mcr 


Mlim 


-(Cp,M)iiiii 


0.45 


2.76639 


0.05 


0.16406 


0.832 


0.890 


0.50 


.50 


2.12.953 


.08 


.27744 


.775. 


.855 


.77 


.55 


1.65519 


.10 . 


.36000 


.742 


.833 


,92 


.60 


1.29190 


.12 


.44856 


.712 


.815 


1.09 


.65 


1.00661 


.15 


.59344 


,.670 


.790 


1,33 


.70 


.77758 


.18 


.75384 


.634 


.760 


1.58 


.75 


.59008 


.20 


.87000 


.610 


.743 


1.74 


.80 


.43381 


.21 . 


.93098 


.598 


.735 


1.82 


.85 


.30.124 


.22 . 


.99396 


.587 


.725. 


1.90 


.90 


.18605 


.25 


1.1953 


.558 


.698 . 


2.15 


.95 


.08783 












1.00 


0 













• TABLE XIII 



VALUES OF THE PRESSURE COEFFICIENT (Cp,M)abs CALCUUTED 
BY imjiS OF EQUATION .(22) 



M 


■ ' 

-(Cp,M)abs 




M . 


1 

-(Cp,M)abs 


0.70 


2.90508 




.. 1.25 


.91103 


.75. 


' . 2.53064 




. 1.30 


.84230 


.80 


2.22420 




. 1.35 


.78106 


.85 


1.97022 




1.40 


.,72627 


.90 


1.75739 




1.45 


. ,67705 


.95 


1,57727 




1.50 


.63266 


1.00 


1.4234-9 




1.55 


.59250 


1.05 


1.29115 




1.60 


,55605 


1.10 


. 1.17644 




1.65 


.52286 


1.15 


■ 1.07636 




1.70 


.49256 


1.20 


.98853' 




1.75 


.46481 



NACA 



Fig.. 1 




Figure 1.- Velocity distribution at the surface of a bump, t = 0.10, for 
several values of the Mach number, 



NACA 



Fig. 2 




Figure 2.- Maximum velocity at the s^irface of a bmip, t = 0.10. 




Figure 3.- Pressure distribution at the surface of a "bump, 
t = 0.10, for M = .83. 



NACA. 



Fig. 




Figare 4.- MaKiraum pressure coefficient as a function of the 
Mach aumter. 



